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High Pressure Insulator-Metal Transition in Molecular Fluid Oxygen 



Marina BastealO, Arthur C. Mitchell and William J. Nellis 
Lawrence Livermore National Laboratory, P. O. Box 808, Livermore, CA 94550 

We report the first experimental evidence for a metallic phase in fluid molecular oxygen. Our 
electrical conductivity measurements of fluid oxygen under dynamic quasi-isentropic compression 
show that a non-metal/metal transition occurs at 3.4 fold compression, 4500 K and 1.2 Mbar. We 
discuss the main features of the electrical conductivity dependence on density and temperature and 
give an interpretation of the nature of the electrical transport mechanisms in fluid oxygen at these 
extreme conditions. 



PACS numbers: 71.30.+h, 62.50.+P, 72.20.-i, 77.22. Ej 
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The transition of condensed matter between electri- 
cally conducting and insulating states is a topic of wide 
scientific interest, whose relevance ranges from supercon- 
ductivity to collosal magnetoresistance JT|, to more re- 
cently thermoelectricity ||. The metal/insulator tran- 
sition has received renewed attention during the past 
decade in the context of high pressure research . Al- 
though much progress has been made in developing ex- 
perimental, theoretical and computational tools appro- 
priate for the study of the metal/insulator transition at 
extreme conditions, our present understanding is mostly 
phenomenological and still incomplete. Given the "sim- 
ple" nature of their interactions, the homonuclear di- 
atomic molecular species, e.g. hydrogen oxygen 
nitrogen |Io) and the halogens pi) , have been in- 
tensively studied using both static and dynamic high 
pressure techniques. 

We report in this letter the first experimental evidence 
for a non-metal/metal transition in the molecular fluid 
phase of oxygen under high dynamic compression. Oxy- 
gen has very rich physics at high pressure. Since dra- 
matic color changes were reported in the lOOkbar pres- 
sure range jl2), the high pressure solid phases of oxy- 
gen have been extensively investigated using structural 
[ fjjl plU , optical |t],H and transport techniques fffl. A 
metallic state and evidence for superconductivity have 
been identified in the solid around IMbar at very low 
temperatures The properties of the liquid on the 

other hand have not been explored much, due to experi- 
mental difficulties and also technical and conceptual chal- 
lenges for theory. 

New insights on the physics of warm dense matter were 
provided by dynamic compression experiments on hy- 
drogen J5jjg]. We present the first electrical conductiv- 
ity measurements of fluid oxygen under dynamic quasi- 
isentropic compression between 0.3 and l.9Mbar. In our 
experiments fluid oxygen reaches up to 4-fold compres- 
sion and temperatures below 7000-fr, conditions never 
before reached experimentally. We note that these con- 
ditions are similar with the ones found in the interiors of 
the giant planets, where oxygen is a major constituent, 
and that these conductivity measurements may be instru- 
mental in explaining the origins of the planetary mag- 
netic fields. 



We measured the electrical resistance of quasi- 
isentropically compressed oxygen starting from high pu- 
rity (99.995%), disk shaped liquid samples [jl6| of density 
do = 1.202(7 / 'cm? at To = 77 K and atmospheric pres- 
sure. The high pressures were generated by multiple re- 
flections of a shock wave between two single-crystal sap- 
phire anvils which enclose the sample. The initial shock 
wave was produced by the impact between a high ve- 
locity projectile onto a stationary target containing the 
sample. The gradual increase in pressure produced by 
this technique yields quasi-isentropic compression. The 
experiments were designed to achieve and maintain ther- 
modynamic steady state conditions at the final pressure 
and density for time durations of the order 100 — 200ns, 
during which the measurements were taken. We opti- 
mized the target geometry and electrodes size and posi- 
tion in order to achieve homogeneous samples of macro- 
scopic size. The shock sensors which provide the trigger 
for the data acquisition were placed outside the sample 
cavity in order to preclude any interference with the mea- 
surements. Other details of the shock reverberation tech- 
nique, cryogenic target construction and electrical cir- 
cuitry were similar to those of Ref. || . 

The equilibrium pressure is determined with a 1% 
accuracy from the measured projectile velocity, using 
the shock impedance matching technique |jl7|| . Densi- 
ties were calculated using a one dimensional hydrody- 
namic code in which the projectile, anvils and sample 
materials were modeled by Mie-Gruneisen and ratio of 
polynomials equations of state (EOS). The EOS's were 
derived based on extensive experimental results of Ref. 
|l8| , ^9| . Temperatures were obtained by chemical equilib- 
rium calculations similar with those in J2(J . The sample 
resistance was extracted from accurate measurements of 
steady voltage and current intensity signals (see Fig.|l| 
- inset, for an example), taken during thermodynamic 
steady state conditions lasting ~ 200ns. The electrical 
resistivity was calculated from the sample resistance us- 
ing R = Cp where C represents the geometrical factor. 
The geometrical factor, or cell constant, was computed as 
a finite element solution to the Maxwell equations for the 
experimental geometry with appropriate boundary con- 
ditions. The calculations spanned a wide range of sample 
thicknesses (between 80(im and 200/im), and were veri- 
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fied by direct measurements performed on materials with 
known conductivities. 

The measured resistivity of fluid oxygen ranges from 
10 9 y£lcm at 0.3A/ bar, down to 8 x W^ttcm a t l.QMbar. 
As seen in Fig.|l| these resistivity values are four to five 
orders of magnitude larger than the resistivities observed 
along the principal shock Hugoniot by Hamilton et al. pl| ] 
at the same pressure. Such a big difference can be ex- 
plained by the fact that the high-pressure states achieved 
using the shock reverberation technique are at signifi- 
cantly lower temperatures (~ 1200X at OAMbar) than 
those on the Hugoniot (~ 6500.fr at OAMbar Q). In the 
single shock (Hugoniot) technique, between 0.18A/6ar 
and OAMbar the densities increase by ~ 20% whereas the 
temperatures nearly triple. It can be inferred therefore 
that the resistivity measured along the principal Hugo- 
niot is mainly due to a much higher degree of thermal 
excitation of the carriers. 

We distinguish two main regimes for the pressure and 
density dependence of the electrical resistivity. At pres- 
sures up to ~ IMbar, there is a rapid drop of the re- 
sistivity of the compressed material which encompasses 
six orders of magnitude. Between ~ IMbar and 2Mbar 
the electrical resistivity shows little sensitivity to pres- 
sure and density variations. In order to understand this 
behavior we need to analyze the changes produced by 
compression on the electronic configuration and chem- 
ical bonding in the warm fluid. The highest occupied 
electronic state in the oxygen molecule at ambient con- 
ditions is ks 12eV below the continuum. In a crystalline 
solid the energy levels of the isolated molecule split and 
widen to form sharply separated, alternating bands of 
allowed and prohibited electronic states. In a fluid, due 
to the inherent structural "disorder" , the energy levels 
broaden to form bands and the band edges spread into 
"tails" . 
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FIG. 1. Pressure dependence of the electrical conductiv- 
ity of shock-compressed fluid oxygen; open circles - new 
quasi-isentropic compression data; filled circles - Hugoniot 
measurements of Hamilton et al.. Error bars are shown where 
larger than symbol size. Inset - example of the measured volt- 
age (lower thick line) and current intensity (upper thin line) 
through the oxygen sample and a parallel 1Q shunt resistance. 



Nonetheless, as our experiments also show, liquid oxy- 
gen at 77 K and lbar is a wide gap electrical insulator. As 
pressure increases the energy bands broaden further and 
create a non-zero electronic density at the Fermi level. 
Most of the states in the tails of the valence and con- 
duction bands and at the Fermi level are localized, as 
demonstrated by Mott p2fl . We postulate the existence 
of an activation energy E a (mobility gap), i.e. the min- 
imum energy that a valence electron needs to acquire 
to move into a delocalized state and participate in con- 
duction, and use an activation model for the electrical 
conductivity variation with temperature T: 
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where a Q and d m are the electrical conductivity and den- 
sity at the transition into the metallic state respectively, 
and Ti.(d m — d) is the Heaviside step function. We as- 
sume a linear dependence of the activation energy on 
the distance between neighboring molecules, Eqn. |2j, and 
zero value above the non-metal/metal transition. This 
approximation should work best in the vicinity of the 
transition, as it represents the first order expansion of 
the mobility gap in the deviation from the electronic 
wave-function overlap at metallization. The fit results 
are shown as triangles in Fig.||. The parameters at 
which the activation energy vanishes are P m — l.2Mbar, 
d m = 'lAg/cm 3 and a = 1205/ilcm. Using the fitted 
a a we extract the activation energy E a — 2fcsTin(cr /cr) 
from the data, Fig.||. Although the activation energy 
effectively vanishes around 1.2Mbar, it becomes in fact 
comparable to the thermal excitation energy at some- 
what lower pressures which is consistent with the ob- 
served change in slope of the experimental data. The 
conductivity at the transition, a , is interpreted as the 
Ioffe-Regel limit |23| , corresponding to a scattering length 
(mean free path) I of the electrons of the order of the in- 
termolecular distance. Using the Drude model, we com- 
pute the effective electronic density n e = a m e v e /e 2 l, 
where v e = hk/m e is the typical electronic speed, kl ~ 1. 
We find that the effective number of conduction electrons 
per molecule is ~ 0.1 (for a good metal, e.g. copper at 
ambient conditions, this number is ~ 0.5). 

The fact that the electrical conductivity does not ap- 
pear to saturate in the metallic state is not unexpected. 
The "resistivity saturation" concept has been introduced 
by Fisk and Webb in order to explain the behavior of 
strongly interacting, disordered systems |24|. However, 
there is both experimental and theoretical |25| evidence 
in support of the absence of saturation. Just as found by 
Millis et al. in the study of a simple representative 
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isentropic compression is mostly molecular. 
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FIG. 2. Density dependence of the electrical conductiv- 
ity of shock compressed oxygen - open and filled circles for 
the present shock-reverberation and Hugoniot experiments of 
Hamilton et al., respectively. Triangles represent the model 
fit (see text) 

system of electrons strongly coupled to phonons [p5| , 
in the case of compressed fluid oxygen there is a smooth 
crossover, as the density is increased, from activated con- 
duction to a transport regime which is dominated by 
hopping between nearest neighbors whose separation dis- 
tance decreases with further increments of pressure. 

In the lower density regime, see Fig.||, there are signif- 
icant deviations of the model from the measured conduc- 
tivities (see also Figj^, the activation energy dependence 
on density) . Linear extrapolation of the oxygen melting 
curve from Ref. Q indicates that up to approximately 
0.55Mbar the states achieved by quasi-isentropic com- 
pression fall into the solid region of the phase diagram. 
We believe that, due to the short time-scales of the exper- 
iments, the system does not crystallize but rather reaches 
a glassy state. The temperature and density calculations 
that we perform using the EOS of the liquid should still 
describe the system reasonably well. However, as the 
oxygen molecule is rather elongated, the glassy states 
so obtained may contain a significant degree of orienta- 
tional order "frozen in" . We speculate that these corre- 
lations lead to increased derealization of the electrons 
and therefore better conduction. 

The model deviations from the single-shock Hugoniot 
conductivity measurements are consistent with our ex- 
pectations that as the density is decreased higher or- 
der terms in the activation energy (E a ) expansion be- 
come important. Also due to the fact that the Hugo- 
niot states are much hotter, it is expected that the fluid 
would dissociate and contributions to the conductivity 
would arise from the atomic species as well. Chemical 
equilibrium calculations similar with those in |20[ yield 
an estimate of the degree of dissociation in the regime of 
interest. For example, in our experiments, at a density 
of 3.9g/cm 3 , IMbar and 3900A approximately 2.4% of 
the molecules break-up. Also, at OAMbar and 1200A" 
the shock-reverberated states have a dissociation frac- 
tion of less than 10 _4 % as compared with the 7% along 
the principal Hugoniot at the same pressure and 6500 A. 
We conclude therefore that the fluid undergoing quasi- 
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FIG. 3. Activation energy (see text) vs. density - present 
data (fluid - open circles, amorphous solid - crossed circles) 
and single-shock data of Hamilton et al. (filled circles) . Errors 
bars are smaller than symbol size except where shown. 

We analyze the crossover to the metallic regime from 
the perspective of the Mott principle, which states that 
an insulator/metal transition may occur when the sep- 
aration distance between the centers, of number density 
n, that produce charge carriers becomes comparable with 
the spatial extent of the electronic wave functions p2| . 
Upon sufficient overlap of the wave functions the out- 
ermost electrons become delocalized and participate in 
conduction. This is predicted to occur at a critical num- 
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a B is a generalized Bohr radius, i.e. the radius where the 
probability of finding a valence electron is maximum. In 
order to estimate a* B for the O2 molecule we calculated 
its electronic charge density as a function of position us- 
ing the GAMESS package. We performed restricted open 
shell Hartree-Fock calculations using Gaussian basis sets 
to determine the electronic structure of the S = 1 ground 
state oxygen molecule. The O2 molecule is 30% elongated 
in the direction of the bond. Since in the warm fluid the 
angular correlations between molecules should be rela- 
tively small, we calculated an effective electronic charge 
density as an average of r 2 ipi(j* over the solid angle. We 
found that for oxygen < r 2 ifjip* > as a function of dis- 
tance has a maximum at 1.2bohr, which was taken to be 
the generalized Bohr radius, a* B . 

The electrical conductivity dependence of liquid oxy- 
gen on the Mott scaling parameter n 1 / 3 x a* B is shown 
in Fig.[|. Also shown are the results for hydrogen, ce- 
sium and rubidium, which are believed to undergo insu- 
lator/metal Mott transitions at high pressures P,^6fl. It 
is remarkable that for the two molecular fluids as well 
as the two alkalis there are obvious changes in the con- 
ductivity at values of the Mott scaling parameter which 
are close to the predicted value for the Mott transition. 
The slower increase of the conductivity in the case of 
the alkali metals, which have only one relatively weakly 
bound valence electron, can be attributed to the spatially 
extended nature of their electronic wave functions. 
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The transition of fluid oxygen to a metallic state should 
be accompanied by the divergence of the dielectric sus- 
ceptibility e, also known as the Goldhammer-Herzfeld cri- 
terion [gTj . Given the limitations of the Clausius-Mossoti 
equation pq], e = (1 + 87rna/3)/(l — Anna/3) (n is O2 
number density and a is the molecular polarizability) , 
we use it only to gain a qualitative understanding. We 
find that at metalization, at a density d = A.lg/cm 3 , 
a w 3.17A 3 , a factor of two bigger than the measured 
polarizability of molecular oxygen at standard ambient 
conditions p9fl . A similar effect has been observed in 
diatomic liquid mercury p(| . We interpret the enhance- 
ment of the O2 molecular polarizability as an indication 
that the molecules at these conditions are distorted. 
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FIG. 4. Electrical conductivity dependence on the Mott 
scaling parameter for oxygen (circles), hydrogen (triangles), 
rubidium (crosses) and cesium (x) . Dotted lines are guides to 
the eye. 

Although, to our knowledge, there are no theoretical 
studies of fluid oxygen in the Mbar pressure range, sev- 
eral predictions have been made p5[ for the solid in con- 
junction with the reported metallization of solid oxygen 
under static high pressure |^||. The reported metalliza- 
tion pressure for the low temperature solid, 0.96Mbar, is 
very close to the onset of the metallic conduction in the 
fluid around \.2Mbar. However, the mechanisms that 
lead to the non- metal/metal transition in the fluid are 
fundamentally different from the solid. In the fluid we 
observe a Mott transition where density driven band clo- 
sure and disorder play the fundamental roles. In the 
solid, on the other hand, the metallization is associated 
with a structural phase transition Jl4| . 
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